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We study in the Ginzburg-Landau approximation, the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) 
phase of QCD with three flavors and one plane wave, including terms of order 0{l/fi). We show that 
the LOFF window is slightly enlarged, and actually splits into two different regions, one characterized 
by u — s and d — u pairings and the other with d — u pairs only. 
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' At high densities and small temperatures quarks are expected to attract each other in the color antisymmetric 
channel. Cooper pairs arc expected to be formed and color superconductivity to arise, see and 01 for reviews of 
' the subject. For three quark flavors, the Color-Flavor-Locking (CFL) phase is known to be the ground state, provided 
\ the density is asymptotically high. The condensate is a spinless color- and flavor- antisymmetric diquark 

■ Pre-asymptotic densities and T ~ arc probably relevant to describe the cores of compact stars (we will consider 
^ I only the zero temperature case in this paper). At those densities one cannot neglect the mass of the strange quark 

■ and the chemical potential differences 5^ due to (3 equilibrium. Much effort has been devoted to study phases which 
might be relevant at such pre-asymptotic densities, such as the 2SC phase 0|, the gapless phases g2SC 0, and, 

, most remarkably, the gapless gCFL phase |3 • The gapless gCFL phase was considered for some time as the most 
' suitable candidate. It was however realized that imaginary gluon Meissner masses (for g2SC see @, and in particular 
\^ \ for the important candidate gCFL see 0) originate instalailities in the gapless phases (and also the 2SC phase shows 
Q • instability H [13 [ill ) . The possibility to remove the instability from the g2SC phase has been recently discussed in 
Ref. [13. 

At all events, for non-zero differences of the chemical potentials, quark pairings with total non- vanishing momenta 
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I. INTRODUCTION 



I ■ may take place, see |l3j and for a review [T3 |. These pairings are called LOFF pairings, from old superconductivity 
Q-i. studies by Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) [13. LOFF pairing gives rise to LOFF phases. In the case of 
^ ' two flavors it has been shown that instability in 2SC indicates that the LOFF phase is energetically favored 



provided there is chromomagnetic stability in the LOFF phase [itI [Tsl . 

The interesting LOFF case for physics is however that of three flavors. In preliminary studies for three flavors it 
was found, in the Ginzburg-Landau (GL) approximation, that condensation of the u — s and d — u pairs is possible in 
the form of inhomogeneous LOFF pairing [l9j, and subsequently it was found that such a phase is chromomagnetic 
' stable [23. The validity of the G-L approximation for such a crystal structure has recently been tested [23|. The 
particular LOFF phase studied in these papers suggested therefore that LOFF phases, in general, could remedy at 
the "impasse" originated from the discovered chromomagnetic instability of gCFL. These studies where made for the 
leading terms in the inverse of the baryon chemical potential and the problem remained to check the validity of such 
approximation. At the same time the color chemical potentials had been neglected on the basis of previous experience 
on the subject, but without quantitative check. In the present paper we complete such investigations by going beyond 
these approximations. 

Recently an important study for three flavors has been completed by Rajagopal and Sharma [23.l23j. always within 
the G-L approximation, but including the study of structures higher than the simplest crystal structure studied by 
us. The authors come out with two preferred structures with face-centered cubic symmetry, one with separate cubic 
structure of u — s and d — u, the other with combined cubic structure for both. Besides the G-L approximation, 
they make use, in their very complex study, of some approximations that we try to overcome in the present paper. 
Although we consider a simpler crystal structure, our results could be relevant also for more complex crystalline 
structures. The other independent aspect, i.e. the validity of the G-L approximation, still remains to be solved, and it 
will probably require long efforts to deflnitely clarify the nature of the stable phases of QCD at intermediate densities 
under neutrality conditions and non negligible strange quark mass. But the importance of comparing with the LOFF 
phases seems unavoidable at this stage of the knowledge, and LOFF phases remain as very serious candidates for the 
solution of QCD at non-asymptotic densities. 
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The plan of the paper is as follows: In section^lwe review the three flavor LOFF phase of QCD. In section UTTl we 
make the expansion. In section HVI we discuss our results. Section IVl contains our conclusions. 



II. REVIEW OF THE THREE FLAVOR LOFF PHASE OF QCD 

In this section we briefly review the main elements for the study of the LOFF phase. The Lagrangean density for 
three flavor ungapped quarks is: 

C = ^P.J^Prf-Mrf + ^^fJo)i'0,■ (1) 



where M"/ = 5"^ diag(0, 0, Ms) is the mass matrix and Dfj^ = dS^^Sij + igAaT^^Sij. The indices i, j refer to flavor 
and a, (3 to color. The matrix /i^^ is a diagonal matrix, which depends on ^ (the average quark chemical potential), 
/ie (the electron chemical potential), and ^3, /is (color chemical potentials) 0. It can be written as follows 
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(i,j = 1,3 flavor indices; a,/3 = 1,3 color indices). Here = idiag(l, — 1, 0), Tg = ^^diag(l, 1, — 2) in color space 

and Q = diag(2/3, — 1/3, — 1/3) in flavor space; /ie is the electron chemical potential; fJ.3,fJ.s are the color chemical 
potentials associated respectively to the color charges T3 and Tg ; /i is the baryonic chemical potential which we fix to 
500 MeV. One has the following chemical potentials for the nine different quarks: 

- _ ^ _ 2ms 
Miib M 2 3 ' 

, Me 2/i8 

Md, = M + ^ , 

, Me 2/i8 

M.S. = M + y - — ■ (3) 

We make use of the High Density Effective Theory (HDET), see HHHIlil and for reviews 0. The quark 
momentum is written as the sum of a large component /in, where n is a unit vector and a residual small component 
i. We also introduce n-dependent fields i/'n and by 

V'(^)=/^e'''— (Vn(^) + *n(x)) ; (4) 

here ipn and correspond respectively to positive and negative energy Dirac solutions. 

We change the spinor basis by defining ipA = Upur,'4'dg,tpbs,'4'dr,^ug,^sr,i'ub,i'sg,ipdb)- This is done through 
unitary matrices Fa, which are reported in Ref. |9|. Finally a Nambu-Jona Lasinio four fermion coupling is added 
to the Lagrangean. This term is taken in the mean field approximation. The procedure corresponds to the same 
coupling and the same approximation as in Ref. 0I- Next, one introduces the Nambu-Gorkov field 

XA = ^(ri. ) , (5) 
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and the Lagrangean can be written in the form 



[ dn r dEd£, , ( {E-^ + Ila)Sab -^ab{v) , 



with {jjL)A — iftu, P-d, fts, P-d, ftu, P-s, P-u, P-s, P-d)- In the above equation E is the energy, ^ = £ • n is the component of 
the residual momentum along n which satisfies |^| < 6, with 6 an ultraviolet cutoff. 
We make a Fulde-Ferrell ansatz for each inhomogeneous pairing. The ansatz is 
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with 



A/(r) = A/exp(2zqi • r) , 



(8) 



where 2qi is the the Cooper pair momentum. In the gap matrix in ® there are three independent gap functions 
Ai(r), A2(r), Aafr). They correspond respectively to d — s. u — s and u — d pairing. The matrix I^ab is given 
expHcitly in 0, 01 ■ 

To determine the vectors qi one has to minimize the free energy. The norms |qi| can be determined by the 
minimization procedure. As for the directions, one separately analyzes different structures to find out which is the 
favorite one. In our previous paper we had neglected the color potentials since they vanish in the normal phase 
and neglected the 0(l//i) corrections. The favored structure was found to be that with the vectors q2 and q3 parallel. 
In the subsequent work |20j it was found that such a phase is chromomagnctic stable. A step forward came from the 
work of Mannarelli, Rajagopal and Sharma who went beyond the G-L approximation confirming its approximate 
validity and the parallel situation of the two vectors in the favorite phase. Here we shall include the 0(1//^) corrections 
and study the role of the color chemical potentials. 

III. EXPANSION 

We assume from the very beginning Ai = and q2 = qs (the hat denotes a unit vector). However we here include 
the chemical potentials /ia and /is. We consider therefore the Ginzburg-Landau functional 



f7 = 17o + f A^ + |a^ + f A^ + |a^ + §^A^A^ + O (A^) . 
where the coefficients aj, /3j and P23 are defined in the following way: 
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/?23 is related to the 3x3 block of the gap matrix, therefore it depends only on Sjiud — Sfiu^^dg and Sfius — SfJ,Ur,sbJ 
the chemical potentials differences are given by 
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We expand the free energy in the chemical potentials {fi^, fig,, fie) around the starting point fi3 — fig — 0, fie 
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where = (^e,^3,/i8) and S^i = ( - ^'/^3,M8 )• Moreover 
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a—r.g.b \ 

is the free energy of the normal phase and 
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Besides the inclusion of color chemical potentials, other effects that might arise from O (l/fi) terms are as follows. 
First we have higher order terms in the expansion of the strange quark energy, which produces the following change 
of the strange Fermi momentum, with respect to the massless case: 
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An effect of the next-to- leading order term in the 1 / expansion is to change the result = — - , which holds in 

. 

the normal phase in the leading approximation. It is worth mentioning here that these corrections do not affect the 
relation between qj and S^j: qj = 1.19975/x/, because this relation is corrected only by terms of the order that 
are negligible within our approximation. 

There are other higher order effects that can be neglected, due to the assumption of weak coupling. First of all we 
can neglect modes outside a shell around the Fermi surfaces, and therefore we assume an ultraviolet cutoff S <^ fi. 
Furthermore we can neglect the contribution of antiquarks, because their effect would produce a correction of the 
order of gfiS where g is the four-fermion Nambu-Jona Lasinio coupling (of dimension [A/]^^). In the weak coupling 
regime g^6 <C g^J,^ ^ 1 and therefore the above mentioned approximation is justified. 



IV. RESULTS 



The neutrality condition 
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with fl expressed by (|18(l determines /ig , /is , /^Sj for which we get the following results. For /is we get 
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We have not included a term proportional to AjA^f^/ie because, as it will be clear below, Sf-ie contains terms of the 
order of at least l//i or A^. For /ig we get 
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Also in this case we do not include terms proportional to A^AgfJ/ig- Finally (S/tg is given by 
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In Eqs. (|23Il and q is defined by q = 1.1997 y/4. On the basis of the results in Eqs. |(52J), we conclude that 
the corrections due to finite chemical potentials and /is correspond to terms proportional to or smaller in the 
free energy. Since in the G-L expansion we take terms up to order 0{A^), we conclude that the chemical potentials 
can be put equal to zero. 

dn dil 

To the equations implementing color and electric neutrality we have to add the gap equations ttt— = 777— = 0. 

In this way we can find the values of A2 A3 , 6^2 , <5/i3 as functions of //i for fixed values of fi { ^ = 500 MeV) and 
the CFL gap Aq (25 MeV in our numerical evaluation). 

The effect of the O (l//i) corrections can be grasped looking at their impact on the distances between the different 
quark chemical potentials. Since fj,^ and /xg are ineffective we consider 
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FIG. 1: Quark chemical potential differences (5/i2 (dashed, blue online) and 5/i3 (dash-dotted, red online) together with their 
ml 

common value — obtained neglecting corrections (full line). 

We note that the effect of introducing the 1 //i correction is to induce an asymmetric splitting of the Fermi surfaces 
of the s and d with respect to the u quark Fermi sphere, more exactly 6^,2 > Sfi^; neglecting higher order effects, 

777, 

as in [13, one would get the result Sfi2 = Sfi^ since in this limit /ig = — -. The results for the splitting of chemical 

4/i 

potentials are reported in Fig. ^ (dash-dotted and dashed, respectively red and blue online) together with the result 
obtained in the large fi limit (continuous black line). 

As a consequence of these results we expect a difference between the two gap parameters, with A2 < A3. This is 
confirmed by the results of the gap equations that are reported in Fig. |21 Also in this case we report the two gaps in 
the new approximation (dash-dotted and dashed, respectively red and blue online) together with the result A2 = A3 
obtained in the large /x limit (continuous black line). 
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FIG. 2: The full line represents the previous solution A2 = A3; the dashed (blue online) and dash-dotted (red online) lines 
are respectively A2 and A3 in the new approximation. All the gap parameters are normalized to the CFL homogeneous gap 
Ao = 25MeV. 




ml/[i (MeV) 



FIG. 3: The free energy difference in the previous and in the present approximation (respectively dashed black line and 
continuous line, red online). 



From these results we can compute the difference between the free energies in the LOFF state and in the normal 
phase. This result is reported in Fig. |31 (continuous curve, red online), together with the result of [Tgl obtained 
neglecting corrections (dotted Hne). 

The effect of the corrections considered in the present paper is to change the nature of the phase transitions in the 
LOFF regime. The LOFF has been enlarged, as it can be seen from both figs. |21 and O However, whereas in the 
previous approximation the LOFF window was characterized by one gap, and by condensation in two channels: u~ s 
and u — d, now we have LOFF pairing in both channels only in a tiny region, for smaller valuers of ml/ ^l, while near 
the end point of the superconductive region {ml/ ^ > 132 MeV) only A3 7^ and therefore there is no u — s pairing. 
We can call this new phase L0FF2s. It differs from the the LOFF phase with two flavors u and d for the presence of 
the strange quark, which plays an active role in modifying, by mass effects, the chemical potentials and being decisive 
to implement electric neutrality. The transition between the two LOFF phases is second order. 

We finally note that we are not able to fix the lower end of the LOFF phase with two gaps A2 7^ and A3 7^ 0. As 
a matter of fact this point would be obtained by comparing the LOFF and the gCFL free energies, which will be only 
feasible going beyond the existing approximation, when all the l/fi corrections in the gCFL phase are considered. In 
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any case we expect that the LOFF phase with two gaps, in the one plane wave approximation, is hmited to a narrow 
range of values of rn^/fi, because, in the previous approximation, its lower end was for rii?^/ ^ 128 MeV [l9l | and 
L0FF2S begins at m^/^ - 132 MeV . 



V. CONCLUSIONS 



An interesting point to be investigated is the chromomagnetic stability of the L0FF2s phase. We have applied the 
results of [23| to this phase, characterized, as we have stressed by A2 = and A3 ^ and we have found that it 
is chromomagnctically stable: three gluons, related to the generators of an unbroken SU(2)c remain masslcss while 
all the Meissner masses of the remainin g glu ons (both the longitudinal and transverse ones) are real, a result that is 
also implicit in the conclusions of refs. Incidentally, since there is change of symmetry, this result shows that 

L0FF2s and the LOFF state with two nonvanishing gaps are distinct phases. It would be tempting to extend the 
analysis of stability to the tiny region with m^/ fi near, but smaller of 132 MeV, and A2 7^ 0, A3 7^ 0. This extension 
is complicated by the fact that we do not know exactly where this LOFF phase extends, because a complete 0(l//i) 
calculation of the gCFL has not yet been performed. We plan to come back to this problem in the future. 

In conclusion we have proved the relevance of effects for the LOFF phase with three flavors in the one-plane 
wav e ap proximation. Some of our results can be important also for the more complex studies with several plane waves 
[22, 123 . Since also in this case the inclusion of the 1/ fi contributions breaks the symmetry between d — u and u — s 
pairings, this can provide some hints on the best strategy to follow in order to identify the most favored crystalline 
structures. Moreover we have found that at the order A^ the color chemical potentials do not vanish, which can 
be also important in crystallography of three-flavor quark matter because in the GL expansion one has to take into 
account terms up to this order. 
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